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pL| ABSTRACT 

' We explore the evolution with redshift of the rest-frame colours and space densities of AGN 

^ hosts (relative to normal galaxies) to shed light on the dominant mechanism that triggers 

accretion onto supermassive black holes as a function of cosmic time. Data from serendipi- 
tous wide-area XMM surveys of the SDSS footprint (XMM/SDSS, Needles in the Haystack 
I— ' survey) are combined with Chandra deep observations in the AEGIS, GOODS-North and 

GOODS-South to compile uniformly selected samples of moderate luminosity X-ray AGN 
' (Lx(2 - lOkeV) ^ lO'^i - lO^^ ergs^^) at redshifts 0.1, 0.3 and 0.8. It is found that the 

fraction of AGN hosted by red versus blue galaxies does not change with redshift. Also, the 
X-ray luminosity density associated with either red or blue AGN hosts remains nearly constant 
since z = 0.8. X-ray AGN represent a roughly fixed fraction of the space density of galaxies 
of given optical luminosity at all redshifts probed by our samples. In contrast the fraction of 
X-ray AGN among galaxies of a given stellar mass decreases with decreasing redshift. These 
findings suggest that the same process or combination of processes for fueling supermassive 
black holes are in operation in the last 5 Gyrs of cosmic time. The data are consistent with a 
picture in which the drop of the accretion power during that period (1 dex since z — 0.8) is 
related to the decline of the space density of available AGN hosts, as a result of the evolution 
of the specific star-formation rate of the overall galaxy population. Scenarios which attribute 
the evolution of moderate luminosity AGN since z « 1 to changes in the suppermassive black 
hole accretion mode are not favored by our results. 
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1 INTRODUCTION 

Understanding how galaxies assemble their stellar mass and how 
the supermassive black holes (SMBH) at their centres grow remain 
major challenges in current astrophysical research. Recent obser- 
vational evidence indicates that the two processes are intimately 
related (e.g. [Ferrarese & MerrittllOOO) [Gebhardt et al.|2000| l and 
therefore should be studied in conjunction: SMBHs affect the evo- 
lution of their hosts and vice versa. Indeed, many modern semi- 
analytic models (SAMs) of galaxy formation include prescriptions 
for the growth of SMBHs (e.g. |Wang et aL||2007[ [Monaco et"ar| 
|2007[ [Somerville et al.||2008| >. In most of those simulations it is 
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assumed that SMBHs, when active, regulate the formation of new 
stars through some form of as yet poorly understood feedback pro- 
cess, thereby affecting the evolutionary path of their hosts. Without 
the energy input from Active Galactic Nuclei (AGN), many SAMs 
fail to reproduce fundamental properties of galaxies, such as their 
colour bimodality (e.g. Cattaneo et al. "2007' '2009) and luminosity 
function (e.g. [Benson et al.(2003 Fontanot et al. 2007b I. As a re- 
sult of the development of AGN/galaxy co-evolution models it is 
now possible to have a physical description of the accretion his- 
tory of the Universe, which observationally, shows a broad plateau 
at 2: ~ 1 — 4 followed by a rapid decline by almost 1 dex to the 
present day (e.g. [Aird et al.[20T0| l. 

The mechanisms adopted by SAMs to trigger AGN activity 
include major gaseous mergers (e.g. |Di Matteo et al.|[2005[ |Cen| 
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|201 \\ , stochastic accretion of cold gas via disk instabilities or mi- 
nor interactions ( [Hopkins & Hemquist|2006 and quiescent inflows 
of hot gas from the galaxy halo l |Croton et a" 2006[ l. Hydrodynamic 
simulations also propose that the secular evolution of galaxies pro- 
vides a large enough reservoir of gas through stellar winds, which 
could give episodic AGN activity and substantial black hole growth 
(e.g. |Ciotti & Ostriker|1997l|2007[ l. To our knowledge this process 
has not as yet been implemented in SAMs. The different SMBH fu- 
eling modes above operate by design under different conditions and 
therefore make specific predictions for the properties of the galax- 
ies that host AGN, including their star-formation rate, morphology 
and environment. This offers an observational tool to test different 
SMBH growth models. 

The evidence above has motivated efforts to study AGN 
hosts using data from large multiwavelength programs, such as the 
AEGIS (All-wavelength Extended Groth Strip International Sur- 
vey; [Davis et a l. 2007 ), the COSMOS (Cosmic Evolution Survey; 
[Scoville et'aL| 2007l, the GOODS (Great Observatories Origins 
Survey) and the XBOOTES ( (Renter et al.|2005[ l. These programs 
use primarily X-rays to locate AGN and multiwavelength obser- 
vations to study their hosts. As a result observational constraints to 
different AGN fueling models have become available in the last few 
years, although the results are often contradictory and the interpre- 
tation may vary among groups. Morphological studies for example, 
have identified a large fraction of disks (~ 30 per cent) among X- 
ray selected AGN at 2 ~ 1, suggesting that major mergers caimot 



be the dominant mode of SMBH growth at that epoch ( Pierce et al 



2007|[Gabor et al.|2009[ [Georgakakis et al.|2009[[Cistemas et al 



201 1 1. X-ray AGN at z « 1 are typically found in small and mod- 



erate size groups ( [Georgakakis et aL 2008 ; Silverman et al. 2009bl 
with dark matter halo masses « 10^^ - 10^^ M© jCoil et al.|2009] 
[Hickox et al.|2009| l. This environment is proposed to be the most 
conducive to merger events ( [Hopkins et al.|2008] but see [Allevato[ 
[et al.|20rT) . Galaxy interactions, particularly those which include 
at least one early-type (red) system, are indeed observed to be more 
common in group environments at z ~ 1 ^Lin et al.|2010| l. The star- 
formation rate of X-ray AGN, although still not well constrained, is 
suggested to depend on their accretion luminosity (e.g. [Lutz et aL] 
[20T0| [Shao et al.[[20T0) . Below the break of the X-ray luminos- 
ity function it is claimed that the star-formation is decoupled to 
AGN activity, consistent with the stochastic accretion mode sce- 
nario. Above the knee of the X-ray luminosity function it is argued 
that the accretion luminosity scales with the star-formation rate, in 
agreement with the major-merger formation picture. 

The observational studies above essentially probe active 
SMBHs close to the peak of the accretion power of the Universe, 
2; « 1, and lack the volume to provide meaningful constraints on 
the AGN population at low redshift z < 0.5. This is an important 
gap in AGN evolution studies, as there are suggestions that the 
dominant mode of SMBH accretion changes with redshift, thereby 
leading to the observed rapid decline of the accretion power at 
z < IjHopkins & Hemquist| [2006[ [Hasinger|[2008l [Lutz et al.[ 
[2010[ >.[Fanidakis et al. (20101 for example, propose two main chan- 
nels of SMBH growth, disk instabilities, which dominate at high 
redshift and bright AGN luminosities, and hot gas accretion, which 
becomes important at low redshift and faint luminosities. In the 
models above shifts between different fueling modes with cosmic 
time should imprint detectable changes in the properties of AGN 
hosts from z — Oto z = 1 and beyond. 

The Sloan Digital Sky Survey (SDSS; [Aba"zajian et al.|2009) 
has identified the largest sample of low redshift (z « 0.1) AGN to- 
date, using diagnostic emission line ratios ( [Kauffmann et al.[2004[ l. 



The selection function of that sample however, is very different 
from that of X-ray AGN in deep surveys, rendering the comparison 
difficult. Recent work has demonstrated that there is only partial 
overlap between AGN samples selected using diagnostic emission 
lines and X-rays (e.g. [Yan et al.|20I I[ [Juneau et al.[[2011^ . This 
might be due to obscuration effects, which are likely to be more se- 
vere at X-rays, dilution of the AGN optical emission lines by stel- 
lar continuum or nebular emission lines from HII regions, intrinsic 
scatter in the X-ray to optical luminosity ratio of AGN or the con- 
troversial nature of sources such as LINERs. This class of objects 
are often included in AGN samples selected by diagnostic emission 
line ratios, although their emission lines may not be dominated by 
accretion onto a SMBH (e.g. |Sarzi et al.|2010[ l. The study of AGN 
hosts as a function of time would benefit from a homogeneous se- 
lection of active SMBHs at all redshifts. 

A step in this direction is the extended Chandra Multiwave- 
length Project (ChaMP, Green et al.| 2009| >, which detects serendip- 
itous X-ray sources on archival Chandra images. That survey cov- 
ers an area of 32 deg"^ and therefore probes a large enough volume 
to detect in a systematic way X-ray AGN at low redshift i Haggard] 
[et al.|[20T0| l. More recently [Georgakakis & Nandra[ ( [201 1 11 com- 
bined archival XMM observations with the SDSS (XMM/SDSS 
survey) to compile a sample of X-ray selected AGN at low red- 
shift, 2 ~ 0.1, over an area of 120 deg^, much larger than the 
extended ChaMP. The advantage of serendipitous X-ray surveys 
like the XMM/SDSS and the ChaMP is that the AGN selection 
function at low redshift is almost identical to that of X-ray AGN 
at 2 ~ 1 and beyond in deep surveys, thereby facilitating the com- 
parison by minimising differential selection effects. [Georgakakis &[ 
[Nandra| ( [20I l\ demonstrated the power of this approach. By com- 
paring X-ray AGN in the XMM/SDSS at 2 ^ 0. 1 with those in the 
AEGIS at 2 ~ 1 they found little evolution of the fraction of AGN 
in red/blue hosts. They argued that this is evidence against a change 
of the dominant SMBH fueling mechanism with redshift. 

In this paper we expand on the [Georgakakis & Nandra[j201 1| | 
results and provide new estimates of the volume density evolution 
of X-ray AGN in red and blue hosts by combining three different 
X-ray surveys with median redshifts 2 = 0.1, 0.3 and 0.8. We also 
determine how the fraction of X-ray AGN in the overall galaxy 
population changes with redshift. These results are discussed in the 
context of semi-analytic models for the growth of SMBHs and the 
evolution of the AGN population. Throughout this paper we adopt 
Ho = 100kms~^Mpc"\ Om = 0.3 and SIa = 0.7. Rest frame 
quantities (e.g. absolute magnitudes, luminosities, stellar masses) 
are parametrised hy h = J/o/100. 



2 X-RAY AGN SAMPLES 

Three different datasets are used to explore the evolution from 
2 ~ 0.1 to 2 ~ 1 of the space density of AGN in blue and red 
cloud hosts. The XMM/SDSS serendipitous X-ray survey ( [Geor-[ 
[gakakis & Nandra[|2011| is used to select AGN at 2 « 0.1. The 
Needles in the Haystack Survey (NHS; Georgantopoulos 2005) is 
employed to study X-ray AGN at a median redshift of 2 ~ 0.3. 
At z ^ 0.8, the AEGIS t pavis et al. 20 07), the GOODS-North 
and the GOODS-South fields are combined to explore the prop- 
erties of the X-ray AGN population. The number of X-ray AGN 
with Lx{2 — lOkeV) > lO*'^ ergs~^ in each of the above three 
samples is presented in Table[T] 
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Figure 1. Sensitivity curves for tlie XMM/SDSS survey (black solid line; 
2-lOkeV band), the NHS (blue dotted line; 2-lOkeV band) and the com- 
bined AEGIS, GOODS-North and GOODS-South fields (red dashed line; 
0.5-lOkeVband). 



2.1 The XMM/SDSS survey 

The XMM/SDSS is a serendipitous XMM survey of the SDSS area 
and is designed to study the properties of X-ray AGN at 2: ~ 0.1. 
Details on the reduction of the XMM observations, the detection of 
sources, the estimation of fluxes and hardness ratios, the determina- 
tion of the survey sensitivity and the identification of X-ray sources 
with optical counterparts are described in |Georgakakis & Nandra| 
(|20TTJ. 

The XMM/SDSS survey includes pointings which have clus- 
ters of galaxies as their prime targets. The overdensity of sources in 
those fields may bias estimates of the fraction of red/blue AGN 
hosts or the fraction of AGN relative to galaxies. We therefore 
exclude from the analysis XMM/SDSS survey fields that have 
clusters as their prime targets. These observations are identified 
from the target name keyword of the event files. This reduces the 
XMM/SDSS survey area to 102 deg^. Sources detected in the 2- 
8 keV band are used. This energy interval is less affected by obscu- 
ration biases compared to softer energies, while the XMM's sensi- 
tivity remains high, thereby resulting in a sample that is sufficiently 
large for statistical studies. The hard band sensitivity curve is plot- 
ted in Figure [T| 

The low redshift X-ray subsample of the XMM/SDSS sur- 
vey consists of 235 hard-band (2-8 keV) detections with 0.03 < 
z < 0.2 and r < 17.77 mag after correcting for Galactic extinc- 
tion ( [Schlegel et al.|1998) . The magnitude cut corresponds to the 
limit of the SDSS Main Galaxy Sample (r < 17.77 mag, |Strauss| 
|et al.|[2002l l, which provides the vast majority of redshifts in the 
SDSS. The photometry is from the New York University Value- 
Added Galaxy Catalog (NYU- VAGC, [Blanton et al.|2005^ which 
uses data from the SDSS DR7 ( Abazajia n et al.|2009 V The NYU- 
VAGC provides better photometric calibration of the SDSS data 
( [Padmanabhan et al.|2008^ compared to the DR7 and an accurate 
description of the SDSS window function. 

XMM prime targets are included in the XMM/SDSS source 
catalogue. They are identified by comparing the positions of the 



detected sources with either the NED coordinates corresponding 
to the target name keyword of the XMM observations or the right 
ascension and declination of the nominal pointing of the XMM ob- 
servations. A total of 42 sources are identified as prime targets in 
the subsample with 0.03 < z < 0.2 and r < 17. 77 mag and are ex- 
cluded from the analysis. We retain however. X-ray sources which 
lie at similar redshifts as the prime target of a particular X-ray ob- 
servation. 

X-ray fluxes are estimated in the standard 2-lOkeV spectral 
band by extrapolating from the observed count rate in the 2-8 keV 
energy interval (see Georgakakis & Nandra 201 1 for details). The 
XMM/SDSS survey detects sources with luminosities as low as 
Lx (2 — 10) ~ 10''" erg s"' (see section 4 for the estimation of X- 
ray luminosity). In contrast, the higher redshift samples (see next 
sections) are sensitive to Lx (2 — 10) ~ 10*' erg s~'. To make fair 
the comparison of AGN at different redshifts we apply a luminosity 
cut of Lx (2 - 10) = 10"' erg s"' and exclude sources fainter that 
this limit. The final XMM/SDSS sample used in this paper consists 
of 175 serendipitous sources with Lx{2 — 10) > 10*' ergs"', 
r < 17.77 mag and 0.03 < z < 0.2 (median redshift z = 0.1). 

2.2 The Needles in the Haystack Survey 

The original NHS is a serendipitous XMM survey which was ini- 
tiated to select and study normal galaxies, i.e. non-AGN, at X-ray 
wavelengths ( [Georgantopoulos et al.|2005| . That survey consisted 
of high Galactic latitude (b > 20) XMM observations that became 
public in June 2004 and overlapped with the second data release of 
the SDSS. The original NHS is now superseded by the XMM/SDSS 
survey (Georgakakis & Nandra|201 l| l. However, extensive optical 
spectroscopy has been obtained for X-ray sources in a subset of the 
original NHS fields (see below). Therefore, we retain that survey 
in its original format but reanalyse the XMM observations using 
the pipeline developed for the reduction of the XMM/SDSS survey 
fields ( [Georgakakis & Nandra|20TT| l. As a result the source detec- 
tion, flux estimation, astrometric corrections, optical identification 
of the X-ray sources and sensitivity map construction in the updated 
NHS follow the methodology described in Georgakakis & Nandra] 
l |20TT) . 

For the subset of the original NHS fields with declination 
S < lOdeg follow-up spectroscopy was carried out at the ESO 
VLT and NTT telescopeaj The total area of that subsample, which 
hereafter will be referred to as NHS, is 5.6 deg^. All X-ray sources 
with extended optical light profiles (i.e. optically resolved) in the 
SDSS and r < 20.5 mag have been observed spectroscopically. 
Because the density of the targets on the sky is low, each source 
was observed separately with a longslit. The observations were car- 
ried out between October 2006 and March 2008 in queue mode 
with the F0RS2 (FOcal Reducer and low dispersion Spectrograph) 
on the VLT (Very Large Telescope) and the EMMI (ESO's Mul- 
timode Instrument) at the NTT (New Technology Telescope). The 
data were reduced using standard routines in IRAF. Redshifts were 
determined for 129 out of 15 1 targeted sources by visual inspection 
of the reduced spectra. The redshift measurements correspond to at 
least two secure spectral features. Figure|2]presents the redshift dis- 
tribution of the sample and demonstrates that at the magnitude limit 
r = 20.5 mag the majority of the sources lie at z < 0.5. 



1 Based on observations made with ESO Telescopes at the La Silla and 
Paranal Observatories under programme IDs 078.B-0623A and 080. B- 
0409A 
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Figure 2. Redshift distribution of NHS X-ray sources with spectroscopic 
follow-up observations from our ESO VLT and NTT programs. 



In the updated NHS catalogue there are 126 hard-band (2- 
8keV) detected sources with r < 20.5 mag and extended optical 
light profiles. Redshift estimates are available for 100 of them (in- 
cluding one Galactic star). The redshift incompleteness is because 
of (i) featureless spectra, (ii) failure to observe some of the targets 
and (iii) mismatches between the updated NHS source catalogue 
and the original one (Georgantopoulos et al. 2005) used to con- 
struct the target list for follow-up spectroscopy. The sample used in 
this paper consists of 73 hard X-ray sources with 0.1 < z < 0.5, 
r < 20.5 mag. The median redshift is z ~ 0.3. The sensitivity 
curve of the NHS in the hard band is plotted in Figure[T] 

Optically unresolved X-ray sources are underrepresented in 
the NHS spectroscopic sample as they have not been targeted by 
our ESO observations. The majority of those sources however, are 
expected to lie at redshift z > 0.5. Out of the 124 unresolved X- 
ray sources with r < 20.5 mag, a total of 28 have redshift avail- 
able from the SDSS. Of them only one lies in the redshift interval 
0.1 < z < 0.5. We therefore expect that the exclusion of optically 
unresolved X-ray sources is not going to have a strong effect on 
the results and conclusions. This issue is further discussed in the 
results section. 



2.3 The deep X-ray surveys 

AGN at z ~ 0.8 are selected from the Chandra surveys of the 
AEGIS, GOODS-North (2 Ms Chandra Deep Field North) and 
GOODS-South (2 Ms Chandra Deep Field South). The Chandra ob- 
servations in the those fields have been reduced and analysed in a 
homogeneous way using the methodology described by |Laird et al.| 
( |2009^ . The optical counterparts of the X-ray sources have been 
identified using the Likelihood Ratio method as described in |Geor-| 
|gakakiseta"L] ( |2009) . For the AEGIS field we use the BRI ground 
based optical photometry obtained as part of the DEEP2 spectro- 
scopic survey ( |Coil et al.|2004| l. For the GOODS North and South 
fields the optical identification of the X-ray sources uses the Hub- 
ble Space Telescope (HST) Advanced Camera for Surveys (ACS) 



imaging observations of those fields (see [Georgakakis et al.|2009| 
for details). 

The spectroscopy for X-ray sources in the AEGIS is from a 
variety of sources: the DEEP2 redshift survey ( [Davis et al.p003l l, 
follow-up observations of X-ray sources at the MMT using the 
Hectospec fiber spectrograph ( [Coil et aL][2009) , the SDSS and a 
number of spectroscopic programs that targeted galaxies in the 
original Groth Strip (Weiner et al.|2005^ . Optical spectroscopy in 
the GOODS-North is available from either programmes that specif- 
ically targeted the X-ray population in these fields (e.g. |Barger[ 
[et al.|2003[|2005l lCowie et al.'2003l) or the Keck Treasury Redshift 
Survey (TKRS;fWirth et al. 2004). There are many spectroscopic 
campaigns in the CDF-South. These include follow-up observa- 
tions targeting specific populations, such as X-ray sour ces (|Szokoly[ 
[et al.|2004l l, JsT-band selecte d galaxies jMig noh et al.|2005 1, and 
high redshift candidates (e.g. [Dickinson et al.|2004(|Stanway et al.[ 
[2004a|b^ , as well as generic spectroscopic surveys of faint galax- 



ies in the GOODS area usmg the FORS2 ( Vanzella et al. 2005 
2008ll and VIMOS ([LeFevre et al.|2004[[Popesso et al.,2008|,2009 



Balestra et al.|2010| l spectrographs at the VLT. 

We select X-ray sources detected in the 0.5-7 keV band with 
optical magnitudes Rab < 24.1 mag for the AEGIS and iab < 
23.5 mag (F775W HST/ ACS band) for the GOODS-North and 
South fields. At this magnitude cut the spectroscopic identifica- 
tion rate of the AEGIS, GOODS-North and GOODS-South X-ray 
sources is 67 (375/557), 85 (118/139) and 87 (105/121) per cent, 
respectively. From the spectroscopically identified sources we se- 
lect those with redshifts in the interval 0.6 < z < 1.2. This results 
in 149 sources in the AEGIS, 66 in the GOODS-North and 54 in 
the GOODS-South. The median redshift of the sample is 0.8. The 
choice of the X-ray band for the detection of sources is to ensure 
that our z ~ 0.8 sample is selected at similar rest-frame energies 
(l-14keV) as the z ^ 0.1 (2-9 keV) and z ^ 0.3 samples (3- 
lOkeV). Differential selection effects are therefore small, thereby 
facilitating the comparison of the AGN properties among the three 
samples. The combined X-ray sensitivity curve for the 3 fields is 
shown in Figure [T] 



3 GALAXY SAMPLES 

In the next sections the space density of X-ray AGN will be 
compared with that of the overall galaxy population. We there- 
fore compile three spectroscopic galaxy samples at z ~ 0.1, 
z ~ 0.3 and z ~ 0.8, which probe similar redshift intervals as 
the low redshift subset of the XMM/SDSS survey, the NHS and the 
AEGIS-l-GOODS X-ray AGN samples respectively. The number of 
galaxies in each sample is listed in Table[T] 

At 2: ^ 0.1, we select a total of 9541 NYU-VAGC sources 
with r < 17.77 mag that overlap with the XMM/SDSS survey 
fields (after excluding galaxy cluster pointings) and have spectro- 
scopic redshifts in the interval 0.03 < z < 0.2. 

The data release 3 of the DEEP2 spectroscopic survey in the 
AEGIS field is used to select galaxies at z fti 0.8. This spectro- 
scopic program uses the DEIMOS spectrograph ( [Faber et al.[2003[ > 
on the 10 m Keck-II telescope to obtain redshifts for galaxies to 
Rab = 24.1 mag. The observational setup uses a moderately high 
resolution grating (R ~ 5000), which provides a velocity accu- 
racy of 30kms~^ and a wavelength coverage of 6500-9100 A. 
This spectral window allows the identification of the strong [O II] 
doublet 3727 A emission line to 2; < 1.4. We use DEEP2 galaxies 
with redshift determinations secure at the > 90% confidence level 
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Table 1. X-ray AGN and galaxy samples 





X-ray 


AGN 








sample 


redshift 


number 


red 


blue 


AGN 








hosts 


hosts 


dominated 


XMM/SDSS 


0.03-0.2 


175 


76 


68 


31 


NHS 


0.1-0.5 


73 


24 


27 


22 


AEGIS/GOODS 


0.6-1.2 


269 


124 


111 


34 


AEGIS/GOODS (stellar mass) 


0.6-1.2 


217 


98 


97 


22 


Galaxies 


sample 


redshift 


number 


red 


blue 










hosts 


hosts 




XMM/SDSS 


0.03-0.2 


8643 


5304 


3339 




AEGIS-NHS 


0.2-0.4 


1960 


291 


1669 




AEGIS 


0.6-1.2 


6797 


1383 


5414 




AEGIS (stellar mass) 


0.6-1.2 


3436 


648 


2788 





The columns are: (1): X-ray AGN or galaxy sample; (2): redshift range of the sources in each sample; 
(3) total number of sources in the sample; (4) number of red cloud sources ; (5) number of blue cloud 
sources; (6) AGN with broad lines and/or and dominant nuclear point source in the optical. 



(quality flag Q > 3; Davis et al. 2007). A total of 6797 galaxies 
with redshifts 0.6 < z < 1.2 are selected to determine the space 
density of the galaxy population at a median redshift of z ~ 0.8 
and to compare against that of X-ray AGN. 

Although the AEGIS DEEP2 spectroscopic survey is geared 
toward high redshift sources, it also includes a large number of 
galaxies at z < 0.6 (4601). We select a total of 1960 galaxies in 
the interval 0.2 < z < 0.4 to estimate their luminosity function at 
z ~ 0.3 and compare it with that of X-ray AGN identified in the 
NHS. 



4 REST-FRAME PROPERTIES 

This section describes how the rest-frame colour, stellar mass, X- 
ray luminosity and intrinsic column density of galaxies and X-ray 
AGN are determined. 

The KCORRECT version 4.2 routines developed by |Blanton &| 
|Roweis| p007j l are used to fit spectral models to the optical pho- 
tometry of X-ray sources and galaxies and then estimate rest-frame 
colours and absolute magnitudes in the AB system. For the sake of 
uniformity, in all three samples used in this paper we estimate rest- 
frame magnitudes in the °'^m, "'^t? bands, which are the SDSS u, 
g filters shifted to z = 0.1. The advantage of this choice of band- 
passes is twofold. They are close to observed frame for our z ~ 0.1 
XMM/SDSS sample and also have effective wavelengths that are 
similar to the rest-frame effective wavelengths of the DEEP2 R 
and / filters at 2: = 1 (Blanton 20061. 

For galaxies selected in the SDSS and X-ray AGN in either 
the XMM/SDSS or the NHS the KCORRECT routines are applied 
to the observed ugriz SDSS photometry. The rest-frame '^'^u and 
°'^g-band magnitudes are estimated from the apparent magnitudes 
of the source in the u and g band filters respectively. In the AEGIS 
field the KCORRECT routines are used to fit models to the BRI 
photometry obtained as part of the DEEP2 survey ( |Coil et al.|2004^ 
and then estimate the rest-frame "'^u and ^'^g-hsad magnitudes 
from the observed DEEP2 R and /-band data respectively. For the 
GOODS-North and South fields the HST/ACS photometry in the 
B (F435W), V (F606W), i (F775W) and 2 (F850LP) bandpasses 
is used as input to KCORRECT. The rest-frame magnitudes in the 
^'^u and "'^g bands are determined from the F435W and F775W 



apparent magnitudes respectively. The "'^(u — <;) vs Mo.ig colour- 
magnitude diagrams of AGN and galaxies at z=0.1 and z=0.8 have 
been presented by |Georgakakis & Nandra| j201 1^ . In that study a 
colour cut of "'^(u — g) ~ 1.4 was employed to separate blue 
from red cloud galaxies. This value is also used here to define AGN 
subsamples associated with red/blue hosts. 

For the estimation of stellar mass, M^tar, we adopt the rela- 
tion between B-band mass-to-light ratio and rest-frame B — R 
colour of [Zibetti et al.| ([2009). We choose to use the mass-to- 
light ratio coefficients estimated for the Jonhson-Cousin filterset 
(see Table Bl of Zibetti et al. 2009) instead of the SDSS bands. 
It was found that the latter set of coefficients, when applied to the 
AEGIS/DEEP2 galaxy sample, produced too many massive galax- 
ies at 2 ~ 0.8 compared to the mass function of Borch et al. 
(2006, see results section). The DEEP2 BRI photometric bands 
are better suited for the estimation of the rest-frame U — B colour 
of galaxies at 2 ~ 1, for which (Zibetti et aL] |2009) do not pro- 
vide mass-to-light ratio coefficients. Therefore for the estimation 
of stellar masses for AEGIS galaxies and AGN we use the ugriz 
photometry obtained as part of the deep synoptic Canada-France- 
Hawaii Telescope Legacy Survey (CFHTLS) to determine their 
B — R colours and B-band luminosities. We use the CFHTLS data 
from the T0003 releastj^ which reaches a limiting magnitude of 
lAB ~ 26.5 mag. The overlap between the CFHTLS photometry 
and the DEEP2 spectroscopic survey is 0.3 deg^. The total num- 
ber of spectroscopically identified galaxies and X-ray AGN in the 
overlap region are listed in Table ^ The adopted approach for de- 
termining stellar masses is not as accurate as fitting galaxy tem- 
plates to the observed spectral energy distribution. Nevertheless, as 
it will be shown in the next sections, it provides stellar mass esti- 
mates which are adequate for our purposes. Also, the results and 
conclusions are not sensitive to systematics (i.e. choise of initial 
mass function) which affect absolute stellar mass estimates. This is 
because the stellar mass distribution of AGN hosts is studied rel- 
ative to that of galaxies (see Results section). The calculation of 
the B — R colour and B-band absolute magnitudes of galaxies and 
X-ray AGN is carried out using the KCORRECT routines. 

The intrinsic column density, Nh, of individual X-ray AGN 
is determined from the hardness ratios between the soft (0.5- 
2keV) and the hard (2-7 keV for AEGIS, GOODS, 2-8 keV for 
XMM/SDSS, NHS) X-ray bands assuming an intrinsic power-law 
X-ray spectrum with index F — 1.9 (e.g. |Nandra & Pounds [1994) 1. 
The derived column densities are then used to convert the X-ray 
flux in either the 2-lOkeV (XMM/SDSS, NHS) or the 0.5-lOkeV 
(AEGIS, GOODS) bands to unabsorbed X-ray luminosity in the 
rest-frame 2-lOkeV energy interval. 

The NHS and AEGIS/GOODS X-ray samples used in this pa- 
per are not sensitive to sources fainter than Lx (2 — 10 keV) = 
10*^ erg s~^. To avoid luminosity dependent biases in the compar- 
ison of the properties X-ray AGN host galaxies at different red- 
shifts we choose to use in the analysis only sources with Lx{'2 ~ 
lOkeV) > 10''^ ergs-^ This criterion excludes 18 XMM/SDSS 
sources. In many studies a brighter luminosity cut is adopted, > 
10*^ erg s~^, to avoid potential contamination of AGN samples by 
normal galaxies. As discussed by Georgakakis & Nandra (2011), 
although normal galaxy candidates are present at faint X-ray lu- 
minosities, AGN remain the dominant population. We provide an 
estimate of the level of contamination by normal galaxies by calcu- 
lating the fraction of X-ray sources with X-ray-to-optical flux ra- 



^ http://terapix.iap.fr/article.php?idarticle=556. 
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tio log/x//opt < —2 and soft X-ray spectra, Nh < 10^^ cm~^. 
These two criteria are often adopted in the hterature for selecting 
normal galaxies in X-ray surveys (e.g |Homschemeier et al.|2003| 
[Georgakakis et al.|2006| l. The fraction of X-ray sources fulfilling 
the above criteria are 3 per cent (5/175) in the XMM/SDSS and 1 
per cent (3/269) for the combined AEGIS and CDF fields. There 
are no sources with log /x/ /opt < — 2 in the NHS. Contamina- 
tion by non-AGN is therefore not expected to affect our results and 
conclusions. 



5 ESTIMATION OF THE SPACE DENSITY OF AGN AND 
GALAXIES 

The space density of X-ray AGN and galaxies in bins of stellar 
mass, (mass-function; MF), Mo.ig absolute magnitude (optical lu- 
minosity function; OLF) and 2-lOkeV luminosity (X-ray luminos- 
ity function; XLF) is derived using the standard non-parametric 
1/Vmax method ( )Schmidt|1968j ). In this calculation we take into 
account the X-ray selection function (in the case of the AGN), the 
optical magnitude limit of different samples and the spectroscopic 
incompleteness. The MF, OLF and XLF in logarithmic bins are es- 
timated by the relation 



V 

* max,! 



(1) 



where Y equals logi\fstar for the MF, Mo-ig in the case of the 
OLF and log Lx (2—10 keV) in the XLF calculation. In the equa- 
tion above Wi is the weight applied to each spectroscopically iden- 
tified source i to correct for the spectroscopic incompleteness of 
the different samples (see below). Vm,ax,i is the maximum comov- 
ing volume for which the source i satisfies the sample selection 
criteria, i.e. redshift range, apparent optical magnitude limit and, 
in the case of X-ray sources. X-ray flux limit. For X-ray AGN the 
Vmax,i depends on X-ray luminosity (Lx), absolute optical magni- 
tude (Mo.ig), redshift (z) as well as the overall shape of the optical 
and X-ray Spectral Energy Distributions (SED) 
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where dV/ dz is the volume element per redshift interval dz. The 
integration limits are zl — zl and z2 = min{zopticai , zu), where 
we have defined zl, zu the lower and upper redshift limits of the 
sample and Zopucai is the redshift at which the source will become 
fainter than the survey optical magnitude limit. ^l{Lx,z) is the 
solid angle of the X-ray survey available to a source with luminos- 
ity Lx and column density Nh at a redshift z (corresponding to 
a flux fx in the X-ray area curve). For galaxies the Vmax,i is a 
function of Mo.ig, z and optical SED only and hence equation |2] 
simplifies to 
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where the symbols are the same as in equation|2]and SI is the survey 
solid angle. The uncertainty at a given luminosity or mass bin is 



= E 



Wi 



14, 



(4) 



Sample variance is not included in the error budget. This is because 
we are mainly interested in the space density of AGN relative to 
galaxies selected in the same field. The effects of sample variance 
are therefore minimised in this differential approach. This is not the 



case for the z — 0.3 sample, for which the X-ray AGN and galax- 
ies are selected from different fields. Nevertheless, the next section 
shows that the AGN XLFs and the galaxy MFs/OLFs determined 
from the samples presented in this paper at median redshifts z=0.1, 
0.3 and 0.8, are in good agreement with previous estimates. This 
further suggests that sample variance effects are small. 

The conversion of the absolute to apparent optical magnitude 
in the 1/Vmax calculation uses the optical k-corrections determined 
by the KCORRECT version 4.2 routines ( [Blanton & RoweiS|2007| l. 
The model that best fits the optical photometric data of a source is 
also used to estimate k-corrections for the same source at different 
redshifts. In the case of the XLF, the intrinsic A'^^r of individual 
X-ray sources is taken into account in the 1/Vmax estimation. The 
X-ray k-corrections are calculated by adopting an absorbed power- 
law spectral energy distribution with F = 1.9 and photoelectric 
absorption cross sections as described by [Morrison & McCammon] 
( 1983 1 for solar metallicity. 

For SDSS galaxies and AGN in the XMM/SDSS survey, the 
weight Wi in equation 1, which corrects for spectroscopic in- 
completeness, is the fraction of the SDSS Main Galaxy sample 
spectroscopic targets that have reliable redshifts in a given sector 
(FGOTMAIN parameter in NYU-VAGC). For the NHS, AEGIS and 
GOODS X-ray sources the weight Wi is estimated following the 
methodology of |Lin et al.]p999) and |Willmer et al.|p006i l. For 
each X-ray source i targeted for spectroscopy the number of X- 
ray sources with (i) reliable redshifts in a given redshift interval, 
Zl < z < zh, (ii) successful redshift determinations outside that 
interval, z < zl oi z > zh, (iii) failed redshifts and (iv) unavail- 
able spectroscopic observations, are counted within spheres defined 
in the three dimensional colour-colour-magnitude space. When es- 
timating the probability, P, that a spectroscopically unidentified 
source lies within the redshift interval zl < z < zh it is assumed 
that unsuccessful redshifts lie outside the redshift range of interest, 
z < Zl ov z > Zh ("minimal" model of [Willmer et al.|2006) . The 
weight Wi for each spectroscopic source is the sum of the probabil- 
ities P of all sources within the colour-colour-magnitude sphere. 
For the redshift limits (zl, zh) we adopt the values (0.1, 0.5) for 
NHS, (0.2, 1.4) for AEGIS and (0.6, 1.2) for GOODS, respectively. 
The data spheres are defined by the g — r, r — i colours and the r- 
band magnitude in the NHS, the B — R, R — I colours and the 
7?-band magnitude in the AEGIS and the B — V , V ~ i and i- 
band magnitude in the GOODS. When estimating the OLF of the 
overall galaxy population using data from the AEGIS DEEP2 spec- 
troscopic survey we include in the weights an additional correction 
which accounts for the probability that a galaxy will be placed on 
a DEEP2 mask (see |Willmer et al.|2006) . 



6 RESULTS 

This section presents the XLF, OLF and MF of X-ray AGN, both 
for the full XMM/SDSS, NHS, AEGIS-l-CDFs samples and for sub- 
samples split by the rest-frame colour of the host galaxy. In this ex- 
ercise we identify and exclude from the analysis sources with either 
broad emission lines in their optical spectra and/or a prominent nu- 
clear point source in their optical images (SDSS or HST/ACS sur- 
vey of the AEGIS/GOODS). The optical continua of those sources 
are contaminated by AGN light. This introduces biases in the deter- 
mination of the host galaxy absolute optical magnitude and stellar 
mass. Also, the colours of those sources do not provide information 
on the dominant stellar population of their hosts. We note however, 
that AGN dominated sources are included in the determination of 



The physics ofAGN evolution 7 



X 
0) 
T3 

O 
ft 

-c; 

OjO 

o 



-s- 

M 
O 



z = 0.8 (AEGIS + CDFs) 



■■35- 



z = 0.4 (NHS) 






1 

■ 


' z=0.1 (XMMSDSS) : 


-e- 
_o 


40 
39 




1 


1 ' 1, ' 1 

z = 0.1 (XMMSDSS) 






1 


1 1 ■ -. 1 




38 
37 




■■'0 

1 


f - 

1 , 1,1" 




42 


43 44 






42 


43 


44 45 46 




-21ogh 


(2-10 keV) (ergs-') 






log 




lOkeV) (ergs-') 



Figure 3. The 2-10 keV X-ray luminosity function. The black filled squares 
are the XLF estimates from the XMM/SDSS survey a\ z ai 0.1 (bottom 
panel), the NHS at z ai 0.3 (middle panel) and the AEGIS+GOODS sur- 
veys at 2 ai 0.8 (top panel). The dotted lines corresponds to the Luminosity 
And Density Evolution (LADE) model of Aird et al. (2010) estimated at the 
median redshift of each sample. In all panels the red open circles are AGN 
in red host galaxies with rest-frame colour " '^(u — g) > 1.4. 
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Figure 4. Space density fraction of X-ray AGN in blue/red hosts as a func- 
tion of lookback time (lower-axis) and redshift (upper x-axis). The vertical 
axis is the ratio of the 2-lOkeV X-ray luminosity function (integrated over 
all luminosities above 10''^ ergs~^) of red (filled circles) and blue (open 
circles) AGN hosts divided by the integrated XLF for all AGN. The blue 
filled circles are offset by -0.2 Gyrs for clarity. The errors are Poisson esti- 
mates and are propagated from the uncertainties of the XLF. 



Figure 5. X-ray luminosity density as a function of hard band luminosity. 
The black filled squares are for the XMM/SDSS survey al z ^ 0.1 (bot- 
tom panel), the NHS at 2 ai 0.3 (middle panel) and the AEGIS-l-GOODS 
surveys at z ai 0.8 (top panel). The dotted lines corresponds to the LADE 
model of Aird et al. (2010) estimated at the median redshift of each sam- 
ple. In all panels the red circles are for AGN in red host galaxies (offset 
by 0.1 dex in luminosity for clarity), the filled blue circles are for blue 
hosts (offset by -0. 1 dex in luminosity) and the cyan triangles are for X- 
ray AGN whose optical colours are dominated by the central engine (offset 
by -0.2 dex in luminosity). The errors are Poisson estimates. 



both the total XLF (i.e. not split by colour) and the total X-ray lu- 
minosity density. 



6.1 The X-ray luminosity function 

Figure [3] and Table |2] present the 2-lOkeV XLFs of the 
XMM/SDSS, NHS and AEGIS-^GOODS surveys at median red- 
shifts of 0.1, 0.3 and 0.8 respectively. At all three redshifts, the 
estimated XLFs are in good agreement with the Luminosity And 
Density Evolution (LADE) model of Aird et al. (2010), which is 
also plotted in the figure. Although our X-ray samples are apparent 
magnitude limited, this selection effect is correctly accounted for in 
the 1/Vmax calculation. For the NHS in particular, the exclusion of 
optically unresolved X-ray sources, for which optical spectra have 
not been obtained as part of our ESO follow-up programs (section 
[2](, does not appear to introduce incompleteness in XLF estimation. 
This is consistent with our finding that optically unresolved sources 
are dominated by AGN outside the redshift range 0. 1 -0.5 probed by 
the NHS. 

Next we split the three samples by rest-frame °'^{u — g) colour 
to explore changes with redshift of the space density of X-ray 
AGN in red/blue hosts and get insights into the dominant accretion 
mode onto SMBHs at different epochs. AGN dominated sources 
are excluded from the analysis. The total numbers of red cloud 
f-^{u ~ g) > 1.4], blue cloud [° '^{u ~ g) < 1.4] and AGN- 
dominated sources in each sample are listed in Table [T] The XLFs 
of red cloud AGN are shown in Figure [3] and are listed in Table |2] 
The ratio of the integrated XLF of AGN in red hosts over the total 
integrated XLF is plotted in Figure |4] It demonstrates that within 
the errors there is no strong evidence for evolution of the red host 
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(X = 13.1, 8 d.o.f.). We therefore find only marginal evidence, 
significant at best at the 95 per cent confidence level (2a for a 
normal distribution), for a change with redshift of the X-ray AGN 
fraction in a given optical luminosity bin. This is surprising given 
that the space density of X-ray AGN drops by almost 1 dex since 
z ^ 0.8. 

6.3 The stellar mass functions of AGN and galaxies 

The conclusions are different if the fraction of galaxies that host 
X-ray AGN is plotted against stellar mass. In this case there is ev- 
idence for evolution in the sense that at higher redshifts a larger 
fraction of galaxies host X-ray AGN at a given stellar mass bin. 
Figure|9]and Table|4]present the mass functions of galaxies and X- 
ray AGN at redshifts z ^ 0.1, 0.3 and 0.8. X-ray AGN with optical 
continua dominated by the central engine (i.e. broad emission lines 
and/or a prominent nuclear point source in their optical image) are 
excluded from the MP estimation. Figure |9] also shows the galaxy 
MPs of |Borch et al.|p006) at similar redshifts. There is fair agree- 
ment with our MF estimates suggesting that our simple approach 
for calculating stellar masses is adequate for statistical studies. The 
agreement with the results of |Borch et al.|p006) is also good if 
the MF is split into red and blue galaxies. Figure |10| plots the ratio 
of X-ray AGN MF over that of galaxies. In all three samples the 
fraction of galaxies that host AGN increases with stellar mass. The 
steepest rise is for the low redshift sample whereas in both the NHS 
and the AEGIS-l-CDFs samples the fraction of X-ray AGN appears 
to flatten above about 10^° Mq. Also, the fraction of galaxies that 
host X-ray AGN drops with decreasing redshift from z — 0.8 to 
z = 0.1. Using the statistical test to compare the z — 0.8 with 
z = 0.1 AGN fraction distribution we findx^ = 31.5(9 d.o.f) and 
a probability P — 99.98 per cent (3.7(7 in the case of a normal dis- 
tribution) that the two distributions are different. This is the result 
of the building up of the mass function of galaxies with time and 
the overall decline of the space density of X-ray AGN to low red- 
shift. At stellar masses > lO'^^ 0, although the statistics are poor, 
there is tentative evidence that the fraction of X-ray AGN is not 
a strong function of redshift. This is likely related to the fact that 
the mass function above this mass limit evolves little with redshift 
since z ^ 0.8 (Borch et al.|2006[|Perez-Gonzalez et al.|2008) . Fig- 
ure [TO] also presents the fraction of X-ray AGN in red/blue hosts. 
The overall trends are similar to the full galaxy and AGN popula- 
tions. Using the test the probabilities that the X-ray AGN frac- 
tion as a function of stellar mass changes from z — 0.8 to z = 0.1 
are P = 0.98 (x^ = 19-3, 9 d.o.f.) for the red AGN subsample and 
P = 0.97 (x^ = 18.8, 9 d.o.f.) for the blue AGN subset. 



7 DISCUSSION 

A slew of AGN fueling modes are proposed in the literature includ- 
ing hot gas accretion (Croton et al. 2006), major galaxy mergers 
(e.g. Di Matteo et al. 2005), internal instabilities or minor interac- 
tions (e.g. Hopkins & Hernquist 2006) and accretion of recycled 
gas from evolved stars (e.g. Ciotti & Ostriker 1997, 2007). Semi- 
analytic models use combinations of these physical processes to 
explain the formation of galaxies and the growth of SMBHs at 
their centres (e.g. ISomerville et al.||2008l IFontanot et al.||2007al 



AGN fraction since z ~ 0.8. Similar results are obtained for the 
AGN in blue hosts. 

More relevant in the study of the dominant AGN accretion 
mode as a function of redshift is the fraction of the accretion power 
of the Universe that is associated with red/blue hosts at different 
epochs. Figure[5]plots the hard band X-ray luminosity density as a 
function of luminosity for the three samples used in this paper. At 
each redshift bin we sum up separately the luminosity density of 
AGN in red/blue hosts and X-ray sources with colours dominated 
by the central engine and divide with the total luminosity density. 
The results are plotted as a function of lookback time and redshift in 
Figure|6] Within the uncertainties there is little change with redshift 
of the relative fraction of the accretion power in red/blue AGN hosts 
from 2 — 0.8 to 2: = 0.1. This statement is quantified by fitting a 
straight line to the luminosity density fraction versus redshift data- 
points of Figure |6] For the slope of the line we find best fit values 
and 68 per cent confidence level errors (Ax' = l)of +O.Olt°;i* 
(blue AGN hosts) and -0.05to:i5 (red AGN hosts). Based on the 
uncertainty in the slope of the linear fit we then estimate that at the 
68 per cent confidence level the luminosity density associated with 
blue AGN hosts cannot increase or decrease more than a factor of 
about 1.4 between z = 0.8 and z = 0.1. For AGN in red hosts we 
find that their luminosity density cannot increase or decrease more 
than factors of 1.6 and 1.2 respectively, from z — 0.8 to z — 0.1 
(68 per cent confidence level). We also caution that the points for 
red/blue AGN hosts in Figure[6]are lower limits as the typical host 
galaxy colours of broad-line QSOs are still not well constrained. 



6.2 The optical luminosity functions of AGN and galaxies 

We also explore changes with redshift of the fraction of X-ray AGN 
relative to the overall galaxy population. For that we estimate the 
optical luminosity function of galaxies and X-ray AGN in Mo-ig 
absolute magnitude bins. X-ray AGN with optical continua domi- 
nated by the central engine are excluded from the OLF calculation 
as their Afo.ig may not be dominated by stellar light. Including 
those sources in the analysis however, would not change the main 
conclusions. The results are presented in Figure|7]and Table|3] For 
comparison, also plotted in that figure are the galaxy OLFs esti- 
mated in Mo.ig bins at 2 ^ 0.1 and 2 ^ 0.8 by| Blanton| j2006^ 
using SDSS-DR4 data and the first data release of the DEEP2 spec- 
troscopic survey respectively. Our galaxy OLFs are in good agree- 
ment with that study. Figure [5] plots the ratio of the X-ray AGN 
OLF over that of the overall galaxy population. Within the errors 
the fraction of X-ray AGN remains roughly constant with redshift 
in the absolute magnitude range shown in the figure. This result 
holds if we split the galaxy and AGN samples by rest-frame colour. 
The x^ test is employed to explore in a quantitative way differ- 
ences between 2 = 0.1 and 2 — 0.8 in the Mo.ig distribution of 
the fraction of AGN. The x' statistic is 

where fi^^^^ i and fi,z=o.8 is the fraction of AGN in the optical 
luminosity bin i at redshifts 2 = 0.1 and 2 — 0.8 respectively, 
and ai is the combined error which accounts for the uncertainty 
in both /i ^=0.1 and /i,z=o.8. For the total X-ray sample we esti- 
mate x' = 15.6 (8 degrees of freedom) which corresponds to a 
probability P — 0.95 that the two distributions are different. For 
the red and blue AGN subsamples the corresponding probabilities 
are P = 0.83 (x' = 11-7, 8 d.o.f.) and P = 0.89 respectively 



[Fanidakis et al.||2010> . In some of these models the evolution 
of the AGN population is intimately related to changes in the 
dominant accretion mode with time (e.g. [Fanidakis et al.[[2010[ l. 
Such changes are expected to imprint on the properties of AGN 
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Figure 6. The fraction of tlie X-ray luminosity density associated with dif- 
ferent AGN samples is plotted as a function of lookback time (lower x-axis) 
and redshift (upper x-axis). The red open circles are for red AGN hosts, the 
blue filled circles represent AGN in blue hosts and the cyan triangles cor- 
respond to AGN with optical colours dominated by the central engine. The 
vertical errorbars are Poisson estimates propagated from the uncertainties 
in the X-ray luminosity density. The horizontal errors represent the redshift 
interval probed by the different samples. For clarity the filled circles and 
triangles are offset by -0.2 and -1-0.2 Gyrs respectively. 
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Figure 7. Optical luminosity function for galaxies and X-ray AGN in the 
XMM/SDSS, NHS and AEGIS/GOODS samples in Mo.ig magnitude bins. 
The red dashed histogram corresponds to X-ray AGN. The blue solid lines 
are for the overall galaxy population. These should be compared with the 
dotted black histograms which are the OLFs of galaxies estimated by Blan- 
ton et al. (2007) at redshifts intervals similar to those of the XMM/SDSS 
and the AEGIS galaxy samples. For the NHS at 2 0.3 there is no esti- 
mate of the galaxy OLF in Ma.ig bins with which we can compare directly 
our results. 
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Figure 8. Ratio of the optical luminosity function X-ray AGN over that of 
galaxies plotted as a function Mo.ig. In each panel the XMM/SDSS, NHS 
and AEGIS/GOODS samples are plotted separately. The top panel is for 
the total sample, the middle is the fraction of red cloud AGN relative to red 
galaxies. The bottom panel plots blue AGN relative to blue galaxies. For 
comparison, in each panel the red dashed histogram shown in the z = 0.3 
and z = 0.8 plots corresponds to the z ^ 0.1 X-ray fraction estimated from 
the XMM/SDSS survey. The errorbars are Poisson estimates propagated 
from the uncertainties in the optical luminosity functions of X-ray AGN 
and galaxies. 
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Figure 9. Stellar mass function of galaxies and X-ray AGN in the 
XMM/SDSS, NHS and AEGIS/GOODS samples. The red dashed his- 
togram corresponds to X-ray AGN. The blue solid lines are for the overall 
galaxy population. These should be compared with the dotted black curves 
which are the Schechter fits to the mass functions of galaxies at ^ = 0, 
z = 0.3 and z = 0.7 presented by Borch et al. (2006) after correcting for 
the different Hubble constant adopted in their paper. 



hosts, e.g. their morphology, star-formation history, stellar mass 
and small/large scale environment. 

In this paper the rest-frame colours of AGN hosts are used as 
proxy for their integrated star-formation history to explore varia- 
tions of this quantity with redshift from z — 0.8 to z = 0.1. Within 
the uncertainties, we do not find evidence that the relative num- 
ber of red/blue X-ray AGN hosts changes with redshift, in agree- 
ment with recent results from ChaMP ( [Haggard et al.|20l0| >. Also, 
the fraction of the accretion density associated with active SMBHs 
in the red sequence and the blue cloud is nearly constant since 
z = 0.8. These findings suggest that the galaxies that host X-ray 
AGN have similar integrated star-formation histories at all redshifts 
probed in this paper, 2: = 0.1, 0.3 and 0.8. This argues against sys- 
tematic changes of the dominant accretion mode from z = 0.8 to 
z — 0.1. The same process or combination of processes for fueling 
SMBHs are likely in operation from z = 0.8 to the present day. 

There are caveats in this conclusion. A large fraction of the ac- 
cretion power at all redshifts is associated with luminous AGN. In 
many of those sources the radiation from the central engine dom- 
inates the optical continuum. The optical colours of those systems 
do not provide any information on the stellar population of their 
host galaxies. Our analysis is therefore insensitive to systematic 
variations with redshift in the fueling mode of broad line AGN. 
Additionally our conclusions are for moderate luminosity AGN, 
Lx (2 - 10 keV) = lO"^ - lO"" erg s"^ and do not apply to QSOs 
with luminosities Lx > 10** ergs~^, as the volumes of all three 
samples used here are too small to identify large numbers of those 
sources. Scattered AGN light may also contaminate the observed 
colours of galaxies, even in the absence of broad optical emission 
lines or a dominant nuclear point sources. Nevertheless, the work 
of [Pierce et al.|p0I0| l suggests that this effect is small. Reddening 
of the intrinsic galaxy colours by dust jCardamone et al.|20I0^ is 
also expected to affect our conclusions. 
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Figure 10. Ratio of the mass function of X-ray AGN over that of galax- 
ies plotted against stellar mass. In each panel the XMM/SDSS, NHS and 
AEGIS/GOODS samples are plotted separately. The top panel is for the 
total sample. The middle panel is the fraction of red cloud AGN relative 
to red galaxies. The bottom panel plots blue AGN relative to blue galax- 
ies. For comparison, in each panel the red dashed histogram shown in the 
z = 0.3 and z = 0.8 plots corresponds to the z ft: 0.1 X-ray fraction es- 
timated from the XMM/SDSS survey. The errorbars are Poisson estimates 
propagated from the uncertainties in the mass functions of X-ray AGN and 
galaxies. 
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If there is little or no change with redshift in the relative impor- 
tance of SMBH fueling modes, how is the strong evolution of the 
AGN population explained? An important piece of evidence for ad- 
dressing this question is the nearly constant fraction of X-ray AGN 
relative to the OLF of galaxies as a function of redshift. This sug- 
gest that the decline of the accretion power of the Universe since 
2 « 1 is linked, at least to the first approximation, to the evolution 
of the galaxy OLF, which in turn is the result of secular evolution 
processes and the decrease with time of the star-formation rate den- 
sity of the Universe. Our results are therefore consistent with a fixed 
accretion mode (or combination of accretion modes) superimposed 
on an evolving galaxy population. 

The stellar mass function of AGN relative to galaxies provides 
complementary information on the evolution of active SMBHs. 
It is found that X-ray AGN hosts out to z ~ 0.8 span a wide 
range of stellar masses, although their fraction increases with stellar 



are likely to be very faint, close t o or below the luminosity limit 
10« 



Kauffmann & Heckman 



2009 1 



mass in agreement with previous studies ( 


Kauffmann et al. 


,2003^ 


Shi et al. 


2008 


[Silverman et al.||2009a| 


Xue et al. 2010 


Hag- 


gard et al. 


2010 


1. This trend is more pronounced at lower redshift. 



where the distribution of the fraction of galaxies with X-ray AGN 
increases nearly monotonically with stellar mass. At higher red- 
shifts however, the distribution levels off at intermediate ma sses, 
> 10^° M©. Similar trends are reported by 



Shi et al. 



(2008 



and 

|Bundy et al.H2008] l, although their X-ray AGN samples are selected 
at brighter luminosities, Lx > 10'*'^ erg s~^, compared to our limit 
of Lx(2 - lOkeV) > lO"^ ergs"^ It is also found that the frac- 
tion of galaxies hosting an active SMBH decreases with redshift 
at a given stellar mass bin. For the red AGN/galaxy populations 
this decrease is because of the build-up of the red galaxy stellar 
mass function from z — 0.8 to 2 — 0.1 and the overall decline 
of the AGN space density. For AGN with blue colours however, 
the blue galaxy MF hardly evolves with redshift (e.g. [Borch et al.| 
|2006[ >. The decline of the fraction of blue galaxies that host AGN 
at a given stellar mass bin is primarily because the space density of 
blue AGN drops with redshift. 

Kauffma nn & Heckman| j2009^ found evidence for two 
regimes of SMBH growth in the nearby Universe. The first ("feast") 
is associated with actively star-forming galaxies and is character- 
ized by a log-normal distribution of accretion rates peaked at a 
few percent of the Eddington limit. For these galaxies the gas sup- 
ply is abundant to feed both star-formation and substantial SMBH 
growth. The second regime ("famine") is associated with evolved 
galaxies and has a power-law distribution of accretion rates. The 
SMBHs of these systems is likely to be fed by recycled gas from 
the evolved stellar population and their growth rates are small. 
The transition between the two modes occurs where the inverse 
of specific star-formation rate, Mstar/SFR, is comparable to the 
age of the Universe, i.e. for galaxies that experience significant 
growth of their stellar mass today. Independent work has shown 
that at a given stellar mass interval the quantity Matar/SFR de- 
creases with increasing redshift (e.g. [Perez-Gonzalez et al.|2008"| l, 
i.e. galaxies were more active in the past. Therefore at a given stel- 
lar mass bin the fraction of galaxies above the Mstar/SFR limit 
which marks the transition from the "famine" to the "feast" SMBH 
growth regimes drops with redshift. If the X-ray AGN selection is 
dominated by SMBH in the "feast" mode then the observed decline 
of the fraction of galaxies that host X-ray AGN from z — 0.8 to 
2 = 0. 1 at a given stellar mass, could be a direct result of the evolu- 
tion of the Mstar/SFR with redshift, i.e. the decline of the volume 
density of potential AGN hosts. 

It is indeed expected that most sources in our X-ray selected 
sample will be in the "feast" mode. AGN in the "famine" regime 



Lx(2- lOkeV) = lU^ergs 
showed that SMBHs that grow in the "famine" regime live in old 
galaxies with 4000 A break amplitude D„(4000) > 1.8. AGN in 
such galaxies typically have [OIII] 5007 emission line luminosi- 
ties, -L[oiiii ^ 10^ — 10^ 1/0 (see Figure 12 of Kauffmann et al. 
2003). Adopting a ratio between X-ray and [OIII] 5007A luminosi- 
ties log(Lx(2 — 10keV)/I/[oiii]) = 1.66 (Georgantopoulos & 
Akylas 2009), we estimate Lx (2 - 10 keV) < 10''^ - 10''^ erg s"^ 
for the optically selected SDSS AGN in the "famine" regime 
(-D„(4000) > 1.8). We caution that there is a large scatter in 
the Lx /-f'[oiii] ratios and faint or obscured AGN may have much 
lower ratios than what is adopte d here (e.g. [Netzer et al.||2006} 



Georgantopoulos & Akylas 2010 1. The conversion of [OIII] 5007A 



to X-ray luminosity is therefore uncertain and the X-ray lumi- 
nosities above may be overestimated by as much as 1 dex. In any 
case these calculations show that most X-ray AGN with Lx (2 — 
lOkeV) > 10"'^ ergs"^ are in the "feast" stage of SMBH growth. 
This is in agreement with the results of [Silverman et aL] l |2009a| > 
who measured 4000 A break amplitudes Z)„ (4000) < 1.8formost 
X-ray AGN at z ^ 1 with Lx{0-5 ~ lOkeV) > lO"*^ ergs"^ in 
the COSMOS survey. 

A number of studies already point to a link between AGN and 
star-formation out to z « 1. The similarities between the global 
star-formation rate density evolution (e.g. [Hopkins|200 4^ and the 
accretion power history of the Universe (e.g. [Aird et al . 2010) sug- 
gests that the two processes are connected, at least in a statisti- 
cal manner (e.g. [Zheng et al.[[2009) . There is also evidence that 
the mean star-formation rate of moderate luminosity X-ray AGN 
hosts evolves with time in a manner that closely mirrors the over- 
all galaxy population ( [Silverman et al.[[2009a| [Lutz et al.][2010[ 
[Shao et al.[2010[ >. This similarity would imprint in the p-band OLF 
of AGN and galaxies, as found in this paper. In apparent contra- 
diction to those results, [Bundy et al.[ ( [200 8) found evidence for 
an association between AGN activity and the quenching of the 
star-formation in galaxies, i.e. their transformation from blue star- 
forming to evolved red systems. However, this does not necessarily 
oppose studies that claim a link between the formation of young 
stars and the growth of SMBHs, as long as the quenching timescale 
is comparable to or longer than the duration of the star-formation 
burst. 

Our analysis suggests that the evolution of the moderate lu- 
minosity X-ray AGN population is related to (i) the presence of a 
specific star-formation rate of galaxies which marks the transition 
from abundant to little gas supply to the SMBH and (ii) the decline 
with redshift of the mean specific star- formation rate of galaxies. In 
this interpretation it is the decrease at lower redshift of the number 
density of potential AGN hosts that drives the evolution of moder- 
ate luminosity X-ray AGN since z = 0.8. This further underlines 
the importance of constraining galaxy formation scenarios to un- 
derstand the growth of SMBHs. 



8 CONCLUSIONS 

Deep small angular size surveys are combined with wide- 
area/shallow samples to study the evolution of moderate luminosity 
X-ray AGN from 2 = 0.1 to 2 = 0.8. The rest-frame colours of 
AGN hosts are used as proxy for their integrated star-formation his- 
tory. 

It is found that the fraction of the accretion density associated 
with red or blue AGN hosts does not change with redshift within the 



1 2 Georgakakis et al. 



uncertainties. This argues against scenarios in wliicii thie dominant 
mode of SMBH growthi clianges witii redsiiift since z — 0.8. 

Tiiere is evidence that the evolution of the AGN population 
is associated with the decrease with redshift of the space density 
of potential AGN hosts. We find that the fraction of galaxies of a 
given g-band absolute magnitude that host X-ray AGN changes lit- 
tle (95 per cent confidence level) from z = 0.8 to z = 0.1, despite 
the rapid decline of the AGN space density at this redshift inter- 
val (close to 1 dex). This suggests a link between AGN, the secular 
evolution of galaxies and the decline of the global star-formation 
rate to the present day. 

When the fraction of AGN among galaxies is plotted against 
stellar mass a systematic decline with redshift is found (99.98 per 
cent confidence level). This is because of the build-up of the stellar 
mass of galaxies and the decrease of the space density of AGN to 
lower redshift. 

It is argued that these trends are consistent with a picture 
where most X-ray AGN live in systems with abundant gas reser- 
voirs, but as the gas supply, and hence the specific star formation 
of the galaxy population, drops with redshift, the number density 
of AGN also decreases. 
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Table 2. X-ray Luminosity Function and X-ray Luminosity Density 



log Lx{2 - 10 keV) 


4,(ALL) 


4>(red) 


4'(blue) 


Lx 4i(ALL) 


Lx <f>{red) 


Lx <i>iblu&) 


(ergs~^) 


(10-* Mpc-^) 


(10-* Mpc-^) 


(10-*Mpc-'') (lO'"* ergs-^ Mpc-'') 


(10^^ ergs-1 Mpc-^) 


(10^* ergs-l Mpc-^) 


XMM/SDSS 


41-42 


4.3 ± 0.7 


2.7± 0.5 


1.6 ± 0.4 


1.1 ± 0.2 


0.6 ± 0.1 


0.5 ± 0.1 


42-43 


1.5 ± 0.2 


0.4 ±0.1 


0.8 ± 0.2 


4.9 ± 0.7 


1.2 ± 0.3 


2.1 ± 0.6 


43-44 


0.14 ± 0.06 


0.08 ± 0.06 


0.012 ± 0.009 


2.7± 1.0 


1.5 ± 0.9 


0.2 ± 0.1 


NHS 


41-42 


9.6 ± 3.5 


4.6 ± 2.3 


5.0 ± 2.6 


3.3 ± 1.1 


1.9 ± 0.3 


1.6 ± 0.7 


42-43 


1.9 ± 0.4 


0.8 ± 0.2 


0.9 ± 0.3 


6.2 ± 1.0 


2.3 ± 0.1 


2.7 ± 0.7 


43-44 


0.24 ± 0.06 


0.04 ±0.02 


0.07 ± 0.3 


6.2 ± 1.7 


0.6 ± 0.2 


1.8 ± 0.9 


AEGIS/GOODS 


41-42 


21.2 ± 5.7 


9.4 ±4.3 


11.7 ± 3.7 


5.3 ± 1.1 


2.5 ± 0.8 


2.8 ± 0.8 


42-43 


5.4 ± 0.7 


2.9 ± 0.4 


2.4 ± 0.5 


20.0 ± 3.4 


10.6 ± 2.6 


9.0 ± 2.3 


43-44 


1.1 ± 0.2 


0.31 ± 0.07 


0.4 ± 0.1 


31.9 ± 4.6 


7.1 ± 1.9 


9.9 ± 2.6 


44-45 


0.12 ± 0.05 


0.02 ± 0.02 


0.01 ± 0.01 


17.7 ±6.5 


1.6 ± 1.6 


3.0 ± 3.0 



The columns are: (1): X-ray luminosity interval; (2): space density for all X-ray AGN; (3) space density for X-ray AGN in red hosts; (4) space density for X-ray AGN in blue hosts; (5) X-ray luminosity density 
for all X-ray AGN; (6) X-ray luminosity density for X-ray AGN in red hosts; (7) X-ray luminosity density for X-ray AGN in blue hosts; 



Table 3, Optical Luminosity Function of AGN and galaxies 









4>AGN(bliie) 


<t>gal(ALL) 




<l>ga,l(bl'Ue) 


(mag) 


(lO-^Mpc-^) 


(10-*Mpc-^) 


(10-* Mpc-^) 


(10* Mpc-^) 


(10* Mpc-^) 


(10* Mpc-^) 


XMM/SDSS 


-18.25 


1.3 ± 0.7 


0.3 ± 0.3 


1.0 ± 0.6 


175.8 ± 7.2 


77.8 ± 4.0 


98.0 ± 6.0 


-18.75 


1.3 ± 0.6 


1.0 ± 0.6 


0.3 ± 0.2 


142.9 ± 4.6 


70.0 ± 2.8 


72.9 ± 3.6 


-19.25 


2.6 ± 0.6 


1.4 ± 0.6 


1.0 ± 0.4 


97.1 ± 2.7 


60.0 ± 1.7 


47.1 ± 2.0 


-19.75 


2.3 ± 0.6 


0.8 ± 0.3 


1.2 ± 0.4 


61.9 ± 1.5 


33.8 ± 1.1 


28.0 ± 1.1 


-20.25 


2.9 ± 0.7 


2.1 ± 0.6 


0.8 ± 0.2 


32.8 ± 0.8 


19.3 ± 0.6 


13.6 ± 0.6 


-20.75 


1.0 ± 0.2 


0.6 ± 0.2 


0.2 ± 0.07 


12.3 ± 0.4 


7.7 ± 0.3 


4.6 ± 0.5 


-21.25 


0.6 ± 0.2 


0.5 ± 0.2 


0.06 ± 0.03 


2.7 ± 0.2 


1.8 ± 0.1 


0.9 ± 0.1 


-21.75 


0.05 ± 0.02 


0.01 ± 0.01 


0.02 ± 0.02 


0.3 ± 0.05 


0.20 ± 0.04 


0.1 ± 0.03 


-22.25 








0.03 ± 0.02 


0.01 ± 0.01 


0.02 ± 0.01 


-22.75 


0.01 ± 0.01 






0.03 ± 0.01 


0.01 ± 0.01 


0.02 ± 0.01 


NHS 


-18.25 


1.4 ± 1.4 




1.4 ± 1.4 


128.8 ± 8.2 


22.6 ± 3.5 


106.1 ± 7.5 


-18.75 


3.6 ± 3.1 


3.2 ± 3.1 


0.3 ± 0.3 


106.3 ± 7.3 


29.1 ± 3.9 


77.3 ± 6.3 


-19.25 


4.7 ± 3.0 


2.2 ± 2.0 


2.4 ± 2.2 


90.6 ± 6.8 


32.6 ± 4.1 


58.0 ± 5.3 


-19.75 


1.4 ± 0.6 


1.1 ± 0.5 


0.2 ± 0.1 


58.2 ± 5.4 


21.6 ± 3.4 


36.6 ± 4.3 


-20.25 


10.4 ± 6.3 


3.8 ± 2.8 


6.5 ± 4.6 


31.8 ± 4.2 


16.6 ± 3.1 


15.2 ± 2.8 


-20.75 


1.1 ± 0.3 


0.4 ± 0.2 


0.4 ± 0.2 


13.2 ± 3.4 


6.1 ± 1.9 


7.1 ± 2.8 


-21.25 


0.6 ± 0.2 


0.2 ± 0.1 


0.02 ± 0.1 


3.1 ± 1.8 


2.1 ± 1.6 


1.0 ± 1.0 


AEGIS/GOODS 


-18.25 


3.2 ± 1.0 


0.4 ± 0.4 


1.2 ± 0.9 


140.2 ± 7.2 


9.4 ±3.2 


130.8 ± 6.6 


-18.75 


7.8 ± 6.2 


1.0 ± 0.6 


6.6 ± 6.1 


124.9 ± 4.9 


17.1 ±2.7 


107.8 ± 4.0 


-19.25 


2.8 ± 1.6 


0.6 ± 0.6 


2.2 ± 1.4 


106.6 ± 3.5 


26.3 ± 2.3 


80.4 ± 2.6 


-19.75 


9.6 ± 3.0 


3.6 ± 1.4 


5.0 ± 2.6 


83.9 ± 2.5 


26.7 ± 1.8 


57.2 ± 1.8 


-20.25 


20.1 ± 8.7 


14.5 ± 8.5 


4.3 ± 1.8 


62.5 ± 1.8 


22.3 ± 1.3 


40.2 ± 1.3 


-20.75 


9.7 ± 2.4 


2.7 ± 1.0 


6.3 ± 2.1 


42.4 ± 1.4 


18.4 ± 1.0 


24.0 ± 1.0 


-21.25 


5.1 ± 1.1 


1.6 ± 0.4 


2.7 ± 1.0 


18.3 ± 0.9 


8.6 ± 0.6 


9.6 ± 0.6 


-21.75 


1.9 ± 0.8 


0.9 ± 0.5 


0.7 ± 0.6 


5.7 ± 0.5 


2.7 ± 0.3 


3.0 ± 0.3 


-22.25 


0.3 ± 0.1 


0.08 ± 0.06 


0.06 ± 0.04 


0.9 ± 0.2 


0.6 ± 0.1 


0.4 ± 0.1 


-22.75 


0.4 ± 0.3 




0.04 ± 0.04 


0.25 ± 0.09 


0.16 ±0.07 


0.10 ±0.06 


-23.25 


0.06 ± 0.04 






0.03 ± 0.03 


0.03 ± 0.03 





The columns are: (1): Absolute magnitude Mq.i^ at the middle of the bin with size AMq .i^ = 0.5.; (2): space density for all X-ray AGN; (3) space density for X-ray AGN in red hosts; (4) space density 
for X-ray AGN in blue hosts; (5) space density of galaxies; (6) space density of red galaxies; (7) space density of blue galaxies. 
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Table 4. Mass Function of AGN and galaxies 



log Mstar 


'I>agn(.^ll) 




<l> AGN (.blue) 




<t>gal('^<^'i) 






(10-*Mpc-^) 


(10-*Mpc-^) 


(10-" Mpc-^) 


(10* Mpc-^) 


(10* Mpc-^) 


(10* Mpc-^) 


XMM/SDSS 


9.1 


0.5 ± 0.6 




0.6 ± 0.6 


412.4 ± 51.3 


203.5 ± 46.7 


208.9 ± 21.1 


9.3 


1.4 ± 1.2 


0.3 ± 0.3 


1.1 ± 1.1 


242.1 ± 19.4 


110.2 ± 15.9 


132.0 ± 11.1 


9.5 


1.6 ± 1.0 




1.6 ± 1.0 


243.9 ± 14.3 


122.8 ± 11.6 


121.1 ± 8.3 


9.7 


2.1 ± 1.0 


0.8 ± 0.8 


1.3 ± 0.6 


263.9 ± 11.2 


159.4 ±9.6 


104.5 ± 6.7 


9.9 


3.5 ± 1.0 


1.6 ± 0.9 


1.9 ± 0.6 


224.4 ± 8.0 


161.5 ± 7.3 


62.9 ± 3.4 


10.1 


r,A ± 1.4 


2.8 ± 1.2 


2.6 ± 0.8 


162.3 ± 5.3 


126.3 ± 4.9 


36.1 ± 2.0 


10.3 


4.5 ± 1.2 


2.8 ± 1.0 


1.8 ± 0.6 


113.7 ± 3.5 


95.0 ± 3.3 


18.7 ± 1.2 


10.5 


4.7 ± 1.4 


4.3 ± 1.4 


0.4 ± 0.1 


66.7 ± 2.1 


59.4 ± 2.0 


7.3 ± 0.6 


10.7 


1.5 ± 0.5 


1.2 ± 0.5 


0.2 ± 0.1 


34.6 ± 1.2 


31.6 ± 1.1 


3.0 ± 0.4 


10.9 


1.8 ± 1.1 


1.8 ± 1.1 


0.04 ± 0.04 


13.1 ± 0.6 


11.9 ± 0.6 


1.2 ± 0.2 


11.1 


0.8 ± 0.4 


0.8 ± 0.4 




3.9 ± 0.3 


3.6 ± 0.3 


0.5 ± 0.1 


11.3 








0.9 ± 0.2 


0.7 ± 0.1 


0.19 ± 0.08 


11.5 


0.02 ± 0.02 


0.02 ± 0.02 




0.3 ± 0.2 


0.3 ± 0.1 


0.05 ± 0.04 


11.7 








0.03 ± 0.03 




0.03 ± 0.03 


NHS 


9.1 








233.7 ±26.0 


29.1pm9.4 


204.6 ± 24.3 


9.3 


0.6 ± 0.6 




0.6 ± 0.6 


153.8 ± 26.0 


33.5 ± 9.7 


120.3 ± 17.6 


9.5 


1.1 ± 1.0 




1.1 ± 1.0 


178.1 ± 22.9 


86.6 ± 17.1 


91.5 ± 15.3 


9.7 


4.8 ± 3.6 


0.4 ± 0.4 


4.4 ± 3.6 


142.0 ± 19.9 


77.9 ± 15.4 


64.1 ± 12.6 


9.9 


11.1 ± 8.9 


2.4 ± 2.1 


8.8 ± 8.7 


100.7 ± 16.0 


61.0 ± 12.5 


39.6 ± 9.9 


10.1 


7.4 ± 5.1 


6.9 ± 5.1 


0.5 ± 0.3 


126.2 ± 18.2 


73.0 ± 13.6 


53.2 ± 12.0 


10.3 


8.6 ± 7.2 


0.8 ± 0.4 


7.7 ± 7.2 


128.5 ± 20.8 


103.5 ± 16.7 


24.9 ± 12.4 


10.5 


1.8 ± 0.9 


1.2 ± 0.9 


0.6 ± 0.2 


57.0 ± 12.7 


51.6 ± 12.1 


5.4 ± 3.8 


10.7 


15.1 ± 10.1 


14.9 ± 10.1 


0.2 ± 0.1 


64.6 ± 14.3 


64.6 ± 14.3 




10.9 


0.8 ± 0.6 


0.6 ± 0.6 


0.2 ± 0.2 


42.6 ± 13.1 


42.6 ± 13.1 




11.1 








4.2 ± 4.2 


4.2 ± 4.2 




11.3 


0.07 ± 0.07 


0.07 ±0.07 










11.5 














11.7 















AEGIS/GOODS 



9.1 








154.6 ± 11.4 


2.9 ± 1.5 


151.7 ± 11.3 


9.3 


6.7 ± 4.4 




6.7 ± 4.4 


172.3 ± 18.6 


8.3 ± 2.9 


164.0 ± 18.4 


9.5 


1.8 ± 0.9 


0.2 ± 0.2 


1.6 ± 0.9 


131.5 ± 8.7 


9.7 ± 2.8 


121.8 ± 8.3 


9.7 


10.2 ± 4.7 


1.1 ± 0.8 


9.1 ± 4.6 


128.1 ± 9.5 


26.6 ±6.8 


101.4 ±6.6 


9.9 


9.6 ± 3.6 


2.8 ± 1.4 


6.9 ± 3.2 


106.7 ±8.3 


32.2 ± 6.2 


74.5 ± 5.5 


10.1 


25.8 ± 16.4 


3.6 ± 2.6 


22.3 ± 16.2 


85.4 ±6.3 


28.5 ±4.5 


57.0 ±4.6 


10.3 


18.0 ± 6.9 


5.3 ± 1.9 


12.8 ± 6.6 


72.1 ± 5.6 


39.0 ± 4.4 


33.1 ± 3.4 


10.5 


19.8 ± 4.6 


10.6 ± 3.7 


9.1 ± 2.7 


72.7 ± 5.4 


42.0 ± 4.5 


30.7 ± 3.1 


10.7 


31.0 ± 21.2 


30.1 ± 21.2 


0.9 ± 0.3 


47.1 ± 3.7 


31.3 ± 3.2 


15.9 ± 1.8 


10.9 


1.7 ± 0.6 


1.1 ± 1.1 


0.6 ± 0.5 


25.8 ± 2.7 


19.3 ± 2.3 


6.4 ± 1.3 


11.1 


2.6 ± 1.1 


2.4 ± 1.1 


0.2 ± 0.2 


17.1 ± 2.3 


12.9 ± 1.9 


4.2 ± 1.4 


11.3 


0.6 ± 0.3 


0.4 ± 0.3 


0.2 ± 0.2 


7.8 ± 1.4 


6.7± 1.3 


1.1 ± 0.6 


11.5 


1.9 ± 1.3 


1.9 ± 1.3 




1.5 ±0.6 


1.3 ± 0.6 


0.2 ± 0.2 


11.7 








0.8 ± 0.4 


0.6 ± 0.4 


0.2 ± 0.2 



The columns are: (1): Stellar mass at the middle of the bin with size A log M^f^^r = 0.2; (2): space density for all X-ray AGN; (3) space density for X-ray AGN in red hosts; (4) space density for X-ray AGN 
in blue hosts; (5) space density of galaxies; (6) space density of red galaxies; (7) space density of blue galaxies. 
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